Multi-state Epidemic Processes on Complex Networks 



Naoki Masuda 1 and Norio Konno 2 

laboratory for Mathematical Neuroscience, RIKEN Brain Science Institute, 
2-1, Hirosawa, Wako, Saitama 351-0198, Japan 
^Faculty of Engineering, Yokohama National University, 
79-5, Tokiwadai, Hodogaya, Yokohama 240-8501, Japan 

Abstract: Infectious diseases are practically represented by models with multiple 
states and complex transition rules corresponding to, for example, birth, death, infec- 
tion, recovery, disease progression, and quarantine. In addition, networks underlying 
infection events are often much more complex than described by meanfield equations 
or regular lattices. In models with simple transition rules such as the SIS and SIR 
models, heterogeneous contact rates are known to decrease epidemic thresholds. We 
analyze steady states of various multi-state disease propagation models with het- 
erogeneous contact rates. In many models, heterogeneity simply decreases epidemic 
thresholds. However, in models with competing pathogens and mutation, coexis- 
tence of different pathogens for small infection rates requires network-independent 
conditions in addition to heterogeneity in contact rates. Furthermore, models with- 
out spontaneous neighbor-independent state transitions, such as cyclically competing 
species, do not show heterogeneity effects. 

PACS numbers: 89.75.Hc, 87.19.Xx, 87.23. Ge, 89.75.Da 

I. INTRODUCTION 

Many epidemic systems can be represented as a graph where vertices stand for individuals 
and an edge connecting a pair of vertices indicates interaction between individuals. Since late 
1990s, real networks underlying disease transmission have been recognized to be complex 
and not approximated by conventional graphs such as lattices, regular trees, or classical 
random graphs. Not only being complex, many networks have the small-world and scale- 



free properties 



An important property of small-world networks is that the distance 



between a pair of vertices, or the minimum number of edges connecting two vertices, is fairly 



small on average. The scale-free property means that the vertex degree k, or the number of 
contacts with other individuals (vertices) for a given individual, has a power law distribution: 
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Pk oc k 7 j3| . Typically, 7 ranges between 2 and 3 |3j, 34]. In contrast, the degree is uniform 
on, for example, the square lattice (k = 4). 



Extending earlier data on heterogeneous sexual contact rates of humans |4Ul6Ll29l]. current 
real data support that scale-free networks underl y propa gation of sexually transmitted dis- 
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eases such as gonorrhea, genital herpes, HIV 
Theoretically, epidemic thresholds above which epidemic outbreaks or endemic states ensue 
scale as (k) / (k 2 ). Here (k) = YlT=i^Pk is the mean degree, and (A; 2 ) = YlT=i^ 2 Pk is the 
second moment of the degree. Then a heterogeneous population, which has a large (k 2 ) 
in comparison with (k), yields a smaller epidemic threshold and are more likely to elicit 
outbreaks. Particularly, scale-free networks with 7 < 3 have epidemic thresholds equal to 
because {k 2 ) diverges. Then, disease propagation on a global scale can occur for an arbitrary 
small infection rate. This holds for the percolation 

called the SIS model) that describes endemic diseases such as gonorrhea 
the SIR model that describes single-shot outbreak dynamics such as HIV 

Apart from heterogeneous contact rates, realistic epidemic processes are also more com- 
plex than merely described by the percolation, the contact process, or the SIR model. 
Practical epidemic models often have more detailed states corresponding to, for example, 
exposed, quarantine, and mutated. Accordingly, state-transition rules become more com- 
plicated It is convenient to represent the transition rules of an epidemic 
model by a schematic figure like Fig. ^ Solid lines mean state-transition routes whose event 
rates are independent of the states of neighbors (e.g. death, recovery, mutation). Dashed 
lines represent ones whose rates are proportional to the neighbors' states (e.g. reproduction, 
infection). Then, the contact process and the SIR model are represented by Fig. E^A) and 
[D^B), respectively. For example, in the contact process, a susceptible individual (state 0) 
turns infected (state 1) at a rate proportional to the number of infected individuals in the 
neighborhood (= n{). The infection rate is denoted by A, and the spontaneous recovery rate 
is set equal to unity. 

Understanding of multi-state contagion models on complex networks may be helpful in 
combatting real epidemics. In this direction, Liu and coworkers analyzed the SIRS model, 
which has the state of partial immunity R and is schematically shown in Fig.^C) 2^j]. They 



also studied the household model, in which each vertex has a graded state corresponding to 
the number of patients in a household [2^. They showed that the epidemic thresholds are 
proportional to (k) / (k 2 ) also in these models. However, how general the universality of the 
(k) I (fc 2 )-scaling is has not been addressed. 

In this paper, we analyze various multi-state epidemic processes originally proposed for 
lattices or populations with homogeneous contact rates. We concentrate on endemic models 
in which the phase boundary between the coexistence state, where both susceptible and 
infected individuals can survive, and the disease-free state with the susceptible individuals 
only is concerned. Of course, we then disregard the other important class of epidemic models 
including the SIR model, in which possibility of one-shot waves of outbreaks rather than 
coexistence is concerned (for this distinction; see QQ). Instead o f covering both classes, 
we investigate phase boundaries of endemic models with more than two states and possibly 
with more than two phases. We show that the critical infection rates are largely extinct 
as (k 2 ) — > 00. However, there are a couple of important exceptions to this rule: models 
with mutation or cyclic interaction. In Sec. |Hl we start with studying two models that 
follow the (k) I (/c 2 )-scaling. Sections ED HV| El and ED are devoted to the analysis of 
representative models with different characteristics and outcomes. Some of them do not 
follow the (k) I (fc 2 )-scaling. The results are summarized and discussed in Sec. IVHI 



II. MODELS WITH SINGLE NEIGHBOR-DEPENDENT TRANSITIONS 

In the following, we investigate various contagion models with different transition rules. 
Because we can never be exhaustive, our strategy is to study known multiple-state models. 
It will turn out that the extinction of epidemic thresholds as (k 2 ) / (k) goes to infinity is 
not a universal story, as summarized in Sec. IVHI What determines the nature of phase 
boundaries is considered to be the gross organization of the transition pathways, some of 
which are independent of the state of the neighbors (solid lines in Figs, d [21 El El Q) and 
others are neighbor-dependent (dashed lines). For the introductory purpose, we analyze in 
this section two three-state models whose epidemic thresholds are proportional to (k) / (k 2 ). 



A. Two-stage Contact Process 



Model 

The state-transition rules of the two-stage contact process |22j are schematically shown 
in Fig. EfA). It differs from the SIRS model (Fig. [D^C)) in that the transition rate from 
to 1 is proportional to the density of the neighbors with state 2 (not 1) of a state-0 vertex. 
This model corresponds to dynamics with two life stages. The states can be interpreted as 
0: vacant, 1: occupied by young individuals, and 2: occupied by adults. Only adults are 
reproductive and generate offsprings in neighboring vacant sites at a birth rate equal to A. In 
other words, a birth event occurs at a vacant site at a rate proportional to A and the number 
of neighboring adults. Youngs (state 1) spend random time of mean r before becoming adult 
(state 2). They are also subject to random death events at a rate of 5. Adults die at a rate 
of 1, which gives normalization of the entire model. Alternatively, we can interpret the three 
states as 0: vacant, 1: partially occupied, and 2: fully occupied colonies. Then, only fully 
occupied colonies are potent enough to colonize vacant lands. In the disease analogue, is 
susceptible, 1 is infected but in the incubation period (not infectious), and 2 is infected and 
infectious. We note that the model corresponds to the SEIRS model without the recovery 
state R. 

The ordinary contact process (Fig. ^A)) is recovered as r = oo. In this state-l 
vertex immediately turns into state 2. The meanfield approximation and the rigorous results 
on regular lattices conclude that there are two phases: {0} (the steady state that consists 
only of state 0) and {0,1,2} (positive probability of steady coexistence of 0, 1, and 2). 
Naturally, large A or r promotes survival [22 1 . Particularly, in the meanfield approximation, 
these two phases are divided by rX = 5 + r. 

Analysis 

To derive meanfield dynamics for populations with heterogeneous contact rates, let us 
denote by pk the probability that a vertex has degree k. Obviously, Y^T=iPk = 1- We denote 
by Pi t k (i = 1, 2) the probability that a vertex with degree k takes state i. The probability 
that a vertex with degree k takes state is equal to 1 — p\ t k — P2,k- The probability that 
a neighbor located at the end of a randomly chosen edge takes state i is denoted by Oj. 
This edge-conditioned probability, particularly 2 in this model, is used to determine the 
effective infection rate 0, Q . 



The dynamics are given by 



p hk = A(l - pi jfc - p2,k)kQ 2 -{5 + r)p 1)k , 
p\k = rp ljk - P2,k- 



(1) 



We note that the effective birth rate (the first term in the first equation of Eq. is 
proportional to /c0 2 , which is the average number of state-2 vertices in the neighborhood of 
a degree-^ vertex. 

When we choose an arbitrary edge, the probability that a specific vertex is connected to 



le probabi lity tna 



this edge is proportional to its degree Therefore, we obtain 

52 k kp k pi, k (t) 



(k) 



(2) 



The steady state is given by 



Pi,k 
* 

P2,k 



P2,k 



XrQ*k 



6 + r + A(r + 1)9^' 

where * indicates the steady state. Plugging Eq. (JHJ) into Eq. (| 

B * _ 1 \- Xr&*k 2 p k 

" 2 (k) ^ 5 + r + X{r + l)e* 2 k' 



(3) 



leads to 



(4) 



Equation (J3J) is satisfied when 02 = 0, corresponding to the disease-free state {0}. When < 
62 < 1, state 2 survives. Equation (JHl) implies pl k ,P2 k > in this situation. Accordingly, 
02 > is equivalent to the {0, 1, 2} phase. Because the RHS of Eq. (j3J) is smaller than the 
LHS when Q 2 = 1> the condition for the {0, 1, 2} phase is 







(b\ 



XrQ^k pk 



dQ* 2 \ (k) V 6 + r + A ( r + l )®*2 k ' 



> 1. 



This yields 



A > 



(S + r) (k) 
r (k 2 ) ' 



(5) 



(6) 



When p k (x k 7 (7 < 3), the critical infection rate vanishes. This agrees with the results for 
the percolation, the contact process, and the SIR model. 



B. Tuberculosis Model 



Model 




A next example is the tuberculosis model shown in Fig. EfB) |45j. In the tuberculosis, a 
majority of the infected individuals does not become infectious before recovery. Accordingly, 
state corresponds to healthy, 1 to infected but not infectious, and 2 to infected and infec- 
tious. This model, which is an extension of the two-stage contact process, has two important 
features. One is that we have two state-transition routes that depend on neighbors' states. 
Both of them depend on the density of infectious individuals (state 2). Antoher feature 
is that disease progression from state 1 to 2 is possible through two parallel routes. An 
early-stage patient (state 1) can develop the disease on its own at a rate of r. At the same 
time, infectious patients (state 2) in the neighborhood increases the possibility of disease 
progression for a state- 1 individual. 

When r = 0, the |0, 1, 2} phase is possible on regular lattices but impossible in perfectly 
mixed populations [45]. When r > 0, the {0} phase and the {0, 1, 2} phase appear for both 
lattices and meanfield populations, depending on values of other parameters. According 
to the ordinary meanfield analysis, these phases are divided by Ar = (5 + r) , which is the 
condition identical to that for the two-stage contact process. 

Analysis 

In heterogeneous populations, the dynamics become 



Pi,k 



A(l - px,k ~ P2,k)kQ 2 - PPi,kk&2 - (5 + r)pi >k , 



p2,k = PPi,kk&2 + rpx.k - P2,k- 



(7) 



The steady state is given by 




5 + r + (X + p + Xr)Q* 2 k + XpQ 



xe* 2 k 




r + pkQ* 2 



1 



(8) 



and 





The {0, 1, 2} phase appears when Q 2 > 0> or 



A > 



(S + r) (k) 
r(k 2 ) 



(10) 



This is identical to Eq. ©• 



III. MODELS WITH COMPETING PATHOGENS 



In many epidemics, multiple strains or pathogens with different transmissivity, virulence, 
and mobility coexist. They compete with each other by preying on common susceptible 
individuals. On top of that, transitions between infected states with different pathogens or 
stages of the disease can occur owing to mutation and disease progression. Coexistence of 
two distinct diseases in a population also introduces competition between them. 

In terms of schematic diagrams, this means that multiple neighbor- dependent infection 
pathways (dashed lines) emanate from a susceptible state. We analyze three models with 
competiting pathogens and show that epidemic thresholds are not entirely governed by the 
degree-dependent factors in two models. 



Let us consider a minimal model for competing pathogens whose transition rules are 
shown in Fig. E[A). State 1 and state 2 correspond to patients infected with strain 1 and 
strain 2, respectively. The two strains compete with each other by feeding on the common 
resource, namely, susceptible individuals (state 0), at rates f3\ and /?2- Any patient sponta- 
neously recovers at rate 1. Unidirectional and spontaneous mutations occur (1 — > 2) at rate 
r [3]. As the mutation rate increases, strain 2 becomes more viable than strain 1. In the 
poupulation ecology context, the three states can be also interpreted as empty (0), species 
A (1), and species B (2). 

Exploiting the fact that the absence of state- 1 vertices is equivalent to the contact process, 
Schinazi showed for regular lattices that there are three phases {0}, {0,2}, and {0,1,2} 
separated by nontrivial critical lines j^J. Small r and large /3i support the {0, 1,2} state. 
The meanfield analysis predicts that the boundary between {0} and {0, 2} is fti = 1, and 
{0, 1, 2} appears when both (3\ > /3 2 ( r + 1) and j3i > r + 1 are additionally satisfied. 

Analysis 

For heterogeneous populations, the dynamics are given by 



A. A Minimal Model with Pathogen Competition and Mutation 



Model 



Pi,k = Pi (1 - Pi,k - P2,k) kQ 1 - (r + l)pi, fc , 




(11) 



The boundary between {0} and {0, 2} is obtained easily. Application of the contact-process 
result implies that {0,2} emerges when 

ft > (12) 

which extends the condition (3 2 > 1 for homogeneous populations. 

We next analyze the boundary between {0, 2} and {0, 1,2}. With the steady state 




k / PiQ* 



+ 1 + (r + 1) ((3 2 Q* 2 + ^ej) k I p ire * + f 3 2 ( r + 1)0* 



(13) 



we obtain 



and 



9 



I x - [3iQ* 2 k 2 p k 

^ l + (B 1 C,+3o)e*k' 1 j 



2 (£;)(r + l)^l + ( / 3 1 C 1 + / 3 2 )e^ 

In the {0, 1, 2} phase, 0i > and 62 > are simultaneously satisfied. The condition 6 2 > 
is satisfied when 

(r + 1) (ifc) 

A > (Jfe2> • (16) 

Equation (JTfij) extends one of the conditions for the {0, 1, 2} phase in the homogeneous 
population: f3\ > r + 1. Imposing Q\ > leads to (3\ > /3 2 (r + l), a condition independent of 
heterogeneity in contact rates. Because of this condition, (k 2 ) — > 00 means survival of state 
2 but not necessarily survival of state 1. Increasing (k 2 ) / (k) makes the survival of state 
1 and state 2 more or less equally likely because they feed on the susceptibles in the same 
manner. Then, the strengths of 1 and 2 must be balanced so that one does not completely 
devour the other. This constrains the range of the mutation rate in a degree-independent 
manner. 

Figure |U shows numerically obtained steady densities of state 1 (A) and 2 (B) in networks 
with 10000 vertices and (k) = 12. Initially, each vertex takes one of the three states inde- 
pendently with probability 1/3. We set (3% = 0.5 and j3 2 = 0.25, which predicts the phase 
boundary: r = (3±/ [3 2 — 1 = 1. 

Scale-free networks with the degree distribution p k oc k 1 with 7 = 2.5 (thickest solid 
lines), 7 = 3.0 (moderate solid lines), and 7 = 4.0 (thinnest solid lines) are produced using 
a static model jlj]. Even though produced networks can be disconnected in general, more 



than 95 % of the vertices constitute one component in every run of our simulations. The 
results for the random graph, whose pk follows the Poisson distribution, are also shown 
(dashed lines). Regardless of networks, the numerical results support the existence of finite 
thresholds in terms of r. The values of critical r are slightly smaller than r = 1 but do not 
vary so much for different pk, as our theory predicts. 



B. Model of Drug-resistant Diseases 

Model 

In many diseases such as tuberculosis, drug-resistant strains emerge by the mutation of 



a wild strain 



Misuse of antibiotics is a major cause of appearance of such stronger 



strains. A four-state model shown in Fig. B) represents spreads of a drug-resistant strain 
(4^ • The states 0, 1, 2, and 3 represent empty, susceptible, infected with the wild strain, 
and infected with the drug-resistant strain, respectively. The wild strain (state 2) and the 
drug- resist ant strain (state 3) compete with each other to devour the susceptible (state 1). 
Only the drug-resistant strain is supposed to be fatal with death rate 1. When infected by 
the wild strain, which is not lethal, an individual recovers at rate r owing to the drug. The 
mutation rate is denoted by <p. 

The model without state 2 (or fl 2 = 0) is equivalent to the SIRS model in Fig. CfC) with 
5 = 0. On regular lattices, {0, 1, 2, 3} appears when r + is small enough. Otherwise, state 



2 is extinguished 



421 ] . Based on the meanfield analysis, the {0, 1, 2, 3} phase ensues when 



>0, /3 2 >r + 0, -xr>&. ( 17 ) 



If any of the conditions in Eq. (JT7j) is violated and /3 3 > 1, then we have {0, 1, 3}. Otherwise, 
only state 1 survives (the {1} phase). 
Analysis 

In heterogeneous populations, we obtain 



Pi,k 

f>2,k 
P3,k 



= (3i(l- p 1>k - p 2 ,k - P3,k) + rp 2 ,k - f32Pi,kk<S> 2 - P3Pi,kkQ 3 , 
= P2pi,kkQ 2 ~ (r + (p)p2,k, 

= /?3Pl,A^©3 + 4>P2,k - P3,k- 



(18) 



The steady densities are given by 



Pl,k 
P*2,k 

V Pik ) 



( 



£1 

A 



\ 



r + (ft 
ftO^ 

\Jft<^ + ft(r + 0) 0* 3 }k J 



(19) 



where 



a = Pi(r + (ft) + [(ft + + A0) fte; + (ft + 1) ft (r + </>) e£] fc. 

Equation (|T9*|) assures that p* & > 0. In addition, 

ft0 



6* 



-e; 



(20) 



(21) 



3 ft - ft(r + W 

provided that ft/(r + 0) > ft. Under this condition, 9^ > is equivalent to the {0, 1, 2, 3} 



phase. Combining Eqs. ©, (fT9*j) and (f2~Tj) results in 



9* 



1 



E 



ftgi^gfe 

A 



(22) 



Since A is of the form: A 
Accordingly, 0?j > when 



C 3 + C^kQ^i Eq. (j22j) is essentially the same as Eq. 
ft(r + 0) (k) 



ft > 



(k 2 ) 



(23) 



Even for large r + <ft, sufficiently heterogeneous pk with large (k 2 ) / (k) makes Eq. (|23|) valid. 
However, for existence of state 2, ft(r + (ft) > ft must be also satisfied. This condition is 
independent of the degree distribution and identical to one of the conditions for the {0, 1, 2, 3} 
phase in the ordinary meanfield case (see Eq. ()17|) ). As is the case for the previous model 
(Sec. IIII A")) , this condition regulates the rates of mutation and spontaneous recovery so that 
the strengths of the wild strain and the drug-resistant strain are roughly balanced. 
If ft/(r + (ft) < ft, then Q* 2 = 0. In this case, Eq. (0 reduces to 



Plk 



ft 



1 



Plk A + (A + i)ftes* \p 3 e* k 



which yields 



3 = 7^£ 



ftfte^v 



(k) V ft + (ft + i)fte^ 



As a result, {0, 1, 3} and {1} are separated by 

8 - {k) 



(k 2 



(24) 



(25) 



(26) 



In summary, the {0, 1, 2, 3} phase dominates even for infinitesimally small infection rates 
when (k 2 ) = oo and /3 2 /(V + 0) > fts- When the latter condition is violated, {0, 1, 3} appears 
for a tiny infection rate. 



C. Superspreader Model 



Model 

Superspreaders, namely, patients that infect many others in comparison with normal pa- 
tients, are identified in the outbreaks in, for example, gonorrhea [if]], HIV [4, 11,0], and 
SARS, |ll 13- Such heterogeneous infection rates are sometimes ascribed to heterogeneous 



16, 



24 



29, 



31. 



461 ] . However, super- 



social or sexual contact rates as specified by pu 4, 
spreaders exist even in diseases whose associated networks are considered to have relatively 
homogeneous k, such as SARS j^. Then, an alternative way to introduce superspreaders 



is to consider two types of patients with different infection rates. Then, competition occurs 
between superspreaders and normal patients. 

We assume that superspreaders and normally infectious patients are not distinguished by 
the strain types. Accordingly, an infection event originating from a superspreader can cause 
a normal patient and vice versa. There is no mutation, and competition occurs indirectly. 
That is how the superspreader model considered here is essentially different from the models 
in Sees. ITTTAl and ITTTRl 

We analyze a three-state model with different infection rates, which is schematically 
shown in Fig. OfC) 18, 44]. The states 0, 1, and 2 mean susceptible, patients of type 1 with 
infection rate Ai and recovery rate 1, and patients of type 2 with infection rate A2 > Ax 
and recovery rate 5, respectively. State 2 corresponds to the superspreader, and in practical 
appliations, we put Ai <C A2. Upon infection, a state-0 vertex changes its state into 1 (2) 
with probability 1 — p (p). There are two neighbor-dependent infection routes in Fig. ED^C). 
State 1 and state 2 indirectly compete with each other through crosstalk. In other words, 
the density of state 1 and that of state 2 affect the rates of transitions — > 2 and — > 1, 
respectively. According to the meanfield analysis, {0, 1, 2} results when 



Ai<5(l -p) + \ 2 p > 5. 



(27) 



Otherwise, {0} results. Slightly weaker but qualitatively similar results have been proved 



for lattices 
Analysis 

In heterogeneous populations, dynamics are given by 

Pi,k = (1 -p) (1 - Pi,k - P2,k) k(XiQi + A 2 2 ) - Pi,k, 

P2,k = P (1 - pi,k - p2,k) k(X 1 Q 1 + A 2 6 2 ) - 5p 2 , k - (28) 

The steady state is: 

pt,k \ (x 1 e* l + x 2 e* 2 )k I s(i - P ) 



Plk ) s + \ P + 6{i- P )]{x 1 e* 1 + x 2 Q^k 



(29) 



V 



which leads to 



e* = W—E1q*. (30) 



Therefore, 6^ > 0, or equivalently, 0* = Ai6^ + A 2 ©2 > 0, is the condition for the {0, 1, 2} 
phase. With Eq. (|2~5j) . we obtain 

0* _ J_ ( 1 ~P) + A 2P] ®*k 2 p k (31) 

' ( k )T S+[p + 5(l-p)]kQ* ' 

and {0,1,2} emerges when 

A 1 5(l-p) + A 2 p>^. (32) 
\k ) 

Equation (|3*2*|) suggests that an arbitrary small infection rate Ai or A 2 allows the endemic state 
as (k 2 ) / (k) — ► oo. In contrast to the models in Sees. IIII AI and IIIIBl degree-independent 
factor does not play a role in this model. This is due to the absence of mutation from 1 to 
2 or from 2 to 1. 



IV. DOUBLE INFECTION MODEL 



Model 

For a realistic network model of populations with natural birth, we have to take into 
account that birth can occur only at empty sites. Then, a minimal model such as the 
contact process must be extended to a three-state process in which states 0,1, and 2 mean 
empty, susceptible, and infected, respectively. A susceptible gives birth to an offspring in a 
neighboring empty vertex, which is operationally similar to a contagion event. An infected 



preys on a susceptible in its neighborhood. Then there are two neighbor- dependent transition 
rates, namely, — > 1 transition at a rate proportional to ri\ and 1 — > 2 transition at a rate 
proportional to n 2 . A similar situation arises in the context of double infection in which 
states 0, 1, and 2 correspond to susceptible, infected, and infected by another pathogen, 
respectively. A second pathogen (state 2) targets individuals infected by a first pathogen 
(state 1). 



15 



Keeping these interpretations in mind, we analyze the model shown in Fig. |3] 
State 2 may more virulent than state 1, and evolution of strains has been investigated with 
a more complicated version of this model Q|. The model extends the contact process and 
the SIR model, which correspond to fi = and A = 5 = r = 0, respectively. 

This model is theoretically intricate. A subtlety stems from the cascaded neighbor- 



dependent transitions — > 1 and 1 



121 ] . There are three basic phases {0}, {0, 1}, 



and {0, 1, 2}. The meanfield solution predicts that {0} and {0, 1} are divided by A = 5. The 
one-dimensional lattice with r = does not permit {0,1,2} p, [jjj. However, on regular 
lattices of any dimension, introduction of a small recovery rate r > elicits {0, 1, 2} j2l|. It 
is counterintuitive in the sense that r is the rate at which state 2 turns into state 1. Indeed, 
the meanfield theory predicts the {0, 1, 2} phase for sufficiently small r fulfilling 

(33, 

Another paradoxical behavior occurs regarding to population density. For example, an 
increase in A lessens the number of state-1 vertices on lattices because a larger A creates 



due to excess mortality 



115, 



more preys to devour for state- 2 vertices 52| . Similarly, too large \x drives state 2 to perish 



391 ]. which defines another {0, 1} phase that we do not examine 
here. This {0, 1} phase is distinct from the {0, 1} phase revealed by the meanfield equations 
and mathematical analysis on lattices. 

This model is complex also in the sense that the phase diagram in the parameter space 
is even qualitatively unknown for lattices. The meanfield solution indicates that the two 
critical lines approach each other as fi — > oo. On the other hand, the y m ay cross at finite A 
and fi, as supported by the improved pair approximation ansatz 
Analysis 



the y m 

HQ- 



In heterogeneous populations, the dynamics read 

Pi,k = A(l - pi }k - p2,k)k®i - 5p ljk + rp 2jk - fipi jk kQ2, 



P2,k = m,fc&©2 - (r + l)p2,fc- 



(34) 



The steady state is given by 

ptk \ = xeik (r + 1 

P * 2k ) s(r + 1) + [A(r + + n&Si k + x^ete* 2 k 2 y ^ 

which leads to 



(35) 



e M J_ v xeik* Pk r + l\ 

Clearly, (6*562) = (0,0), which corresponds to the {0} phase, solves Eq. (}3T)j) . To explore 
other phases, let us set 



Jl{ 19 2 ' ~ (k) ^ 5{r + l) + [\{r + l)Q\ + ixQl}k + \^Q\Qlk 2 ' K ' 

and 

f (&* - A ^ Q l V ^ 1 (-1Q) 

J2y lJ 2, ~ (k) ^ 5{r + 1) + [A(r + 1)61 + ^©2] k + XpQ\Q* 2 k 2 ' 1 ; 

Then Eq. ® is equivalent to 6*/i(6*, 0*) = and 0^f 2 (0i, ©2) = °- 

First, we identify the boundary between {0} and {0, 1}. In the {0, 1} phase, we have 
©2 = 0. Under this condition, we look for < Q\ < 1 that satisfies /i(9*,0) = 0. By 
substituting 62 = into Eq. (J3*7|) . we obtain 

= ^Er£fe- 1 = a (39) 

k 

Because 

. . A v-^ k 2 pu A k 2 pu , , 

/•M^EjTS-is^ElF- 150 ' (40) 

Eq. ({HOI) is satisfied when /i(0, 0) > 0, that is, 

- > (41) 

This conclusion complies with the results for the contact process (Sec. HJ). 

Second, we examine the boundary between {0, 1} and {0, 1, 2}. The {0, 1, 2} phase implies 
6^ > and 63 > such that /1 (61,62) = and / 2 (@!,©2) = °- Let us suppose that 
Eq. (j4*Tj) is satisfied because {0} results otherwise. Then, there is a < 61 < 1 that satisfies 



Eq. flU. Since df x /dQ\ < and dfx/dQ* < 0, there is a curve /i(9*, 9*) = in the 9*-9* 
space that looks like the solid line (I) (when /i(0, 1) < 0), or (II) (when /i(0, 1) > 0) in 
Fig. HA). 

Because df2/Q\ > and 8/2/0^2 < 0> h monotonically increases in Q\ on lines 9^ = 
and 6* + 9a = 1. If the isocline /2(9*, <d* 2 ) = is located as the dashed line in Fig. E^B), 
it nontrivially crosses fiiQ^Q^) = ( one °f ^ ne solid lines) to yield the {0,1,2} phase. 
Because /2(C), 9£) = — 1 (0 < 9^ < 1), ^2(1, 0) > is necessary for the dashed line to exist. 
Then, there is a unique < 9 < 1 so that / 2 (9, 1 - 9) = 0. Finally, /i(9*, Q* 2 ) = and 
/ 2 (9*,92) = cross when (i) / 2 (9i,0) > and (ii) A(9,l - 9) < are both satisfied. 
Using 

hOdu 0) = ^ ^ - 1 = 0, (42) 



A 



the condition (i) is equivalent to 



r (g o] - XfiQl V k3pk 1 

/2(Ol ' 0) " (r + l)(A;)^^9^" 1 



A^9x ^ k 3 Pk fx5 - 



(r + 1) (k) 5 + XQik (r + l)A* 



^ -fl + 7 ^J >0, (43) 



(r + l)(k) V (r + l)A 
that is, 

A + ^r<W' (44) 

which is an extension of the ordinary meanfield solution (Eq. ([33)1 ). We note that (i) implies 
^(1,0) > 0, which we required beforehand. 
To show (ii), define 



We obtain 



A(r + l)9 + /i(l-9) 
/i(9,l-9) 

= ^(9, 1-9) + 0^(9, 1-9) 

____ _A_ y> (r + 1 + C^Qk)k 2 p k 

(k) 2f S(r + 1) + [A(r + 1)9 + (l - 9)] fc + A//9 (l - 9) k 2 5 

(r + 1 + C 5f iQk)k 2 Pk i _ c 

(k) ^ [A(r + 1)9 + (1 - 9)] k + A/i9 (l - 9) P 5 



To summarize, both {0} and {0, 1} disappear to be replaced by {0, 1, 2} as (k 2 ) — ► oo. 



V. ROCK-SCISSORS-PAPER GAME 



Model 

Winnerless cyclic competition among different phenotvpes abounds in nature. For ex- 
ample, real microbial communities of Escherichia coli |l9i . l2(l ] and color polymorphisms of 



natural lizards 



481 ] have cyclically dominating three states. The evolutionary public-good 
game with volunteering (choice of not joining the game) also defines a three-state popula- 
tion dynamics with cyclic competition |47l. |49j . A typical consequence of these dynamics is 
oscillatory population density with each phenotype alternatively dominant. 

Let us consider a simple rock-scissors-paper game with cyclic competition (Fig.0A)). The 
ordinary meanfield analysis yields a single phase: a neutrally stable periodic orbit in which 
the densities of states 0, 1, and 2 are cyclically dominant 11711 . The coexistence equilibrium 



inside the limit cycle is stabilized on regular lattices [13|, |51( and trees . Consequently, 
convergence to {0, 1, 2} with a damped oscillation occurs on these graphs. We focus on the 
steady states only in the following analysis. 
Analysis 

The dynamics in heterogeneous populations are given by 

Pi,k = A(l - pi jk - P2,k)k®i - Wi,kkQ2, 

p2,k = PPi,kkQ 2 - P2,kk{l - ©i - 2 ), (47) 



which yields 




xei I i-ei-e 



2 



(Ae; + ^es)(i - e$ - e* 2 ) + A/ie^ \ ^ 

Noting that Eq. (jjHj) is independent of k, we obtain 



(48) 





X+p+1 I A 

for any k. The degree distribution does not affect the steady state. 



(49) 



Similarly, it is straightforward to show that degree distributions are irrelevant in the voter 
model (Fig.E^B)) and cyclic interaction models with more than three states 
from these examples is that heterogeneous contact rates do not influence the equilibrium 
population density if there is no neighbor- independent state transition. However, models 
with at least one neighbor-independent transitions show extinction of the epidemic threshold 
as (k 2 ) I (k) — > oo. For example, variants of the rock-scissors-paper game supplied with 
spontaneous death or mutation rates .52] are expected to belong to the class analyzed in 
Sec. EH 




. A lesson 



VI. TWO-POPULATION MODEL 



Model 

In many infectious diseases of humans, such as malaria, yellow feber, and dengue feber, 
transmission are mediated by other hosts such as mosquitoes. In these diseases, direct 
human- human or mosquito-mosquito infection is absent, and two distinct populations trans- 
mit diseases between each other. As an example, we analyze a model for malaria spreads 
(Anderson and May, 1991, Ch. 14). Following the original notation, the meanfield dynamics 
are represented by 

y = ^(1 -y)ip -7y, 

if) = acy(l — if;) — flip, (50) 

where y and if) are the proportions of infected humans and mosquitoes, respectively. The 
size of the human population is denoted by N, and Q is the size of the female mosquito 
population; only female mosquitoes infect humans. In addition, a is the bite rate, b and 
c are infection rates, 7 is the recovery rate of the humans, and /x is the death rate of the 
mosquitoes. Infected humans and infected mosquitoes survive simultaneously if they do. 
Based on Eq. J5DJ), this occurs when 

a 2 bcVt 

B ' S N^ >1 - (51) 

Otherwise, the malaria is eradicated. 
Analysis 



How does this condition change by introducing heterogeneous contact rates? In network 
terminology, relevant networks are bipartite, with each part representing the human pop- 
ulation and the mosquito population. Let us denote the degree distributions of humans 
and mosquitoes by {p y ,k} and {p^ t k}, respectively. Real data do not suggest that {p y ,k} 
or {Pip,k} is scale-free. However, we treat general degree distributions because the obtained 
results may apply to other multi-population contagion processes. 

The dynamics are given by 

Vk = -j—{l-y k )kB^-jy k , 

ip k = ac(l - ip k )kO y - fiipk, (52) 
where y k {ipk) is the probability that humans (mosquitoes) with degree k are infected, and 



y k 



©V = 77T- 5^ ^.fc^fe, ( k )^ = ^2 k P^,k- (53) 
k k 

The steady state is calculated as 



/ abVLQ^k ac&* y k \ 

{VI rk) = [N 1 + abne;k> fi + ac&;k ) ' (54) 

which is solved trivially by (6*, O^,) = (0, 0). To explore nontrivial solutions, let us eliminate 
from Eqs. (JHHJ) and (JSH) to obtain 



The RHS of Eq. (|55|) is less than 1 when 0* = 1. The endemic state results if 



d .„ a 2 bcQ (k 2 ) y (k 2 )^ 



dQ* y 



(RHS of Eq. (jSSD) = AT 1 1 \ ii \ > h (56) 

which extends Eq. (j51|) . Equation (J51)|) indicates that divergence of just either (k 2 ) y or 
(k 2 )^ is sufficient for the epidemic threshold to disappear. Even if both moments are finite, 
their effects are multiplicative. This is consistent with the results for two-sex models of 
heterosexual HIV transmissions j^J and the bond percolation on bipartite graphs [31I l^. 



VII. CONCLUSIONS 



We have analyzed various models of endemic infectious diseases in populations with het- 
erogeneous contact rates, or on complex networks. The effects of heterogeneity on epidemic 
thresholds are summarized in Tab. In many models, diverging second moments of the 
degree distribution extinguishes the epidemic thresholds, as reported previously for the per- 
colation I9 , llil 33], the contact process 16 1 




SIRS model 25], and the household moodel 



ioias, as reported previously ior tne 
M El, the SIR model Q, H Q, the 



261 ] . On scale-free networks, which underly 



sexually transmitted diseases and computer viruses (see Sec. Hlfor references), introduction 
of a tiny amount of virus to a population can cause an endemic state even in models with 
more complex transition rules. 

However, the models with competing pathogens and mutation fSecs. ITTL^I and lTH B|) show 
different behavior. The heterogeneity in contact rates equally boosts the infection strengths 
of competing pathogens in these models. Whether multiple strains survive or one overwhelms 
the others depends on network-independent mutation rates that modulate relative strengths 
of pathogens. The rock-scissors-paper game and the voter model stipulate another class of 
models in which heterogeneity does not affect the equilibrium population density at all. For 
the heterogeneity to take effects, there must be at least one neighbor-independent transition 
rate. 

For a fixed graph or a contact rate distribution, one can create new models that are not 
convered by this paper. Our analysis has not been exhaustive. However, we consider that 
what essentially matters is gross arrangements of contagion pathways. If there is at least one 
transmission route independent of the neighbors' states, which is not the case for the rock- 
srissors-paper game, the equilibrium population density will depend on degree distributions 
as does the contact process. A care must be paid to cases of competing pathogens with 
mutation. We hope that our results give prescription for understanding complex-network 
consequences of other models. 

A last note is on the stability of solutions. For the contact process on complex networks, 
the coexistence phase is stable if it exists j^, [3]. Similarly, the stability of the simplest 
nontrivial phase is assured for models analyzed in this paper. However, stability analysis 
of other phases seems mathematically difficult. In addition, the stability of coexistence 
solutions of the rock-scissors-paper game shows network dependence j^. These topics are 



warranted for future work. 
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FIG. 1: The state-transition rules of (A) the contact process, (B) the SIR model, and (C) the SIRS 
model. Solid lines represent transitions independent of states of neighbors. Dashed lines represent 
neighbor-dependent transitions. The values indicate the transmission rates, where is the number 
of neighbors with state i. 
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FIG. 3: The state-transition rules of the models with competing pathogens analyzed in (A) 
Sec. Em (B) Sec. IfflBl and (C) Sec. lmTl 
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FIG. 4: Numerically obtained stationary densities of (A) state 1 and (B) state 2 for the model in 
Sec. IIII Al (see Fig. EfA) for the transition rules). The scale-free networks with 7 = 2.5 (thickest 
solid lines), 7 = 3.0 (moderate solid lines), 7 = 4.0 (thinnest solid lines), and the random graph 
(dotted lines) are used. See the text for other parameter values. 




FIG. 5: The state-transition rules of the double infection model. 




FIG. 6: (A) The isocline f\ = (solid lines) and sign of fi, (B) f'2 = (dashed line) and sign of /2 
superimposed on f\ = (solid lines). 



TABLE I: Summary of the effects of heterogeneous contact rates on epidemic thresholds. 



model 


phases 


threshold — > 

as {k 2 ) — * oo 


remarks 


section 


CP 


{0}, {0,1} 


yes 


similar for SIR, 
SIRS, household 


in 


two-stage CP 


{0}, {0,1,2} 


yes 




urn 


tuberculosis 


{0}, {0,1,2} 


yes 




HLU 


2 pathogens 
with mutation 


{0}, {0,2}, 
{0,1,2} 


yes 


{0, 1, 2} requires a 
Pfc-indep. condition 


|11IA| 


drug-resistant 


{1}, {0,1,3}, 
{0,1,2,3} 


yes 


{0, 1, 2, 3} requires a 
Pfc-indep. condition 


|111 H| 


superspreader 


{0}, {0,1,2} 


yes 




Unci 


double infection 


{0}, {0,1}, {0,1,2} 


yes for both 
{0} and {0, 1} 




ED 


rock-scissors-paper 


{0,1,2} 


no 


similar for 
the voter model 


El 


malaria 


nonendemic, endemic 


yes 




ED 



